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Abstract 



Radiative-recoil corrections to hyperfine splitting in muonium of orders 
a(Za)(m/M) 2 Ep and (Z 2 a)(Za)(m/M) 2 Ep are calculated. These correc- 
tions are of the second order in the small electron-muon mass ratio. An 
analytic expression [(—6 In 2 — j)a(Za) — j^(Z 2 a)(Za)](jj) 2 Ep is obtained. 
Numerically the correction is equal to —0.0351 kHz and is of the same order 
of magnitude as the expected accuracy of the current Los Alamos experiment 
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I. INTRODUCTION 



The hyperfine splitting in the ground state of muonium has been measured with high 
precision |j] 

Az/ Mu (n = 1) = 4 463 302.88 (16) kHz 5 = 3.6 ■ HT 8 , (1) 

and there are further opportunities for improving the accuracy. The lifetime of the 
higher hyperfine state is extremely large r = 1 ■ 10 13 s and gives negligible contribution to 
the linewidth. Thus the linewidth is completely defined by the muon lifetime r M « 2.2- s 
which leads to the width T^/h = 72.3 kHz. The experimental accuracy attained thus far 
corresponds to measuring the energy splitting with an experimental error of Ai> exp / (T^/h) m 
2.2 • 10~ 3 of the natural linewidth. A new experiment is now running at Los Alamos 0, 
with the goal of measuring Ah>u u { n — 1) with a precision of about 1 • 10~ 8 , or at the level of 
a few hundredths of a kilohertz. 

In anticipation of the long-awaited results from this precision experiment, recent years 
have witnessed a surge of theoretical activities in this field. A tedious calculation of the 
purely radiative corrections of order a 2 (Za)E F [|]-|9|] was followed by works [|TQ|jT| JTT| T p~2 



where leading corrections of higher orders (a(Za) 3 lnZa Ep, a 2 {Za) 2 \v? Za Ep, 
a 2 (Za) 2 \nZa E F , (Za) 3 (m/M) In 2 Za Ep, (Z a) 3 (m / M) \n Za\n(m/M)E F , 
a(Za) 2 (m/M) In 2 Za Ep) were obtained. Some other corrections of higher orders 
(a 2 (Za){m/M) ln 3 (m/M) Ep, a 2 (Za)(m/M)\n 2 (m/M) E F ) have already been known for 
some time [|T3"| , |PI| | . The magnitude of these higher order corrections varies from a few thou- 
sandths to a few tenths of a kilohertz and some of them turned out to be rather large. The 
correction of order a(Za) 2 E F was also recalculated recently p5|JT6|| . 



In this paper we consider radiative-recoil corrections of order a(Za)(m/M) 2 E F and 
(Z 2 a)(Za)(m/M) 2 E F , which were not discussed so far. For muonium both the fine struc- 
ture constant a and the electron-muon mass ratio (m/M) are of the same order of magnitude, 
and one could anticipate that the corrections under discussion would produce contributions 
comparable to the higher order corrections of order a 4 Ep mentioned above. Another rea- 
son for interest in the radiative-recoil correction of second order in mass ratio is connected 
with a long standing discrepancy between the results of |T7[ and ]T8"|,|H| for radiative-recoil 



corrections of the first order in mass ratio, induced by the radiative insertions in the electron 
line. The discrepancy is about three standard deviations of the accuracy of the numerical 
calculations in |I7| |, and is equal to 0.22 kHz. This discrepancy, which was of a purely aca- 



demic interest when the respective calculations were performed, now acquires some practical 
importance. An apparent reason for this contradiction could be connected with the fact that 
while in an explicit expansion in the small mass ratio was performed, with only the 



contribution of the first order preserved, in the calculations in [[HJ only the nonrecoil con- 
tribution was thrown away, and consequently the result of this work contains terms linear 
in the mass ratio as well as all contributions of higher orders in the mass ratio. With the 
results of this work we will be able to test this hypothesis about the origin of the discrepancy 



between the results in |17j and |lq . [19| . 

Radiative-recoil corrections to hyperfine splitting are generated by the diagrams with 
all possible radiative insertions in the two-photon exchange graphs. It suffices to calculate 
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respective matrix elements between the Coulomb-Schrodinger wave functions (see, e.g., dis- 
cussion in lpO|| ). We will consider below separately the corrections generated by the three 
types of diagrams: polarization insertions in the exchanged photons, radiative insertions in 
the electron line, and radiative insertions in the muon line. 



II. POLARIZATION OPERATOR CONTRIBUTIONS 

The general expression for the a[Za)Ep corrections to hyperfine splitting induced by 
the electron vacuum polarization insertions in the exchanged photons has the form (see, e.g., 

) 

= a(Za) E F -L / * / ^jf&r (2) 
7r u jo J m z k A — k [1 — v ) 



\k 2 + fx^ko + iO k 2 + iOJ k 2 -2k 

where the exchanged momentum k is measured in units of the electron mass and the 
dimensionless parameter /i is given by the expression fi = m/ (2M). Note that the expression 
in Eq. (0) differs from the integral for the skeleton graph with two exchanged photons only 
by the substitution 

J_ _ - C A ^ 2 (1 ~^ 2 /3) 

P^Wo V k 2 (l-v 2 ) -4' ^ 

Let us now describe briefly the calculation of the contribution of relative order fi 2 con- 
tained in Eq. (|2|). First we combine the electron and the photon denominators, and obtain 

8E^ = a(Za)E F (--L) C t^~ 2 v 2 (l - ^) C dz (4) 
v Ti /i ' Jo 1 — v v 3 ' Jo 

d , fd 4 k 1 3^-2k 2 



v Oa 2 ' J t7c 2 k*-fi- 2 k 2 { -k 2 + 2b p k + a 2 p Y 
where a 2 = 4z/(l — v 2 ), b p = 1 — z. 

The general expression in Eq. (§]) contains nonrecoil radiative corrections of order 
a(Za)Ep pi,22l, as well as radiative-recoil corrections of all orders in the electron-muon 
mass ratio generated by the polarization operator insertions and admits numerical calcula- 
tion. Due to the structure of the integrand explicit extraction of the corrections of definite 
order in the mass ratio is more involved. Direct application of the standard Feynman pa- 
rameter methods leads to such integrals for the radiative-recoil corrections which do not 
admit expansion of the integrand over the small mass ratio before integration, making the 
analytic calculation virtually impossible. Instead of combining denominators in Eq. (f|) with 
the help of the Feynman parameters we use an approach developed earlier for calculation of 



the radiative-recoil corrections of order a(Za)(m/M)Ep (see, e.g., review in p0[). The idea 
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is to perform integration over the exchanged momentum directly in spherical coordinates. 
Following this route we come to the expression 



SEP* = a{ Za) E F (-^) j 



dv 



IT 2 ' JO 1 — V 2 



1 - 



dz 



(5) 




k 2 (k 2 + a 2 ) 



+ a 



2\2 



2%(k) + ^(k) 



which readily admits expansion in the small parameter \i. We have introduced here stan- 
dard auxiliary functions <& n (k), which emerge in all calculations in this paper. Definitions 
and properties of these functions are described in Appendix [A|. 

The crucial property of the momentum integrand in Eq. (|5|), which facilitates further 
calculation, is that the denominator admits expansion in the small parameter \x prior to 
momentum integration. This is true due to the inequality 4 / u 2 fe 2 ,&; 4 / [k 2 + a 2 ) 2 < 4fi 2 which is 
obviously satisfied because < a p < oo and < b p < 1. In this way we may easily reproduce 



nonrecoil radiative corrections induced by the vacuum polarization insertions |21| , |22| , and 
respective corrections of the first order in the mass ratio which were obtained in [^,^4]]. We 
will omit here these calculations and describe only calculation of the corrections of order 
(m/M) 2 . Expanding the denominator in Eq. (0) in the small mass ratio we obtain 



5E^ « a(Za)E F (-%) C - 

v 7T 'JO 1 



dv 



- V 



v 2 [l 



— / dz 
3 / Jo 



(6) 



2^{k) + 2^(k) + &£(k) + 2^\k) + $f (k) } j. 

Functions &f(k) do not depend on the small parameter \x (see Appendix |A|) and, hence, 
do not generate corrections of order fi 2 . Calculation of the /^-integrals (integrals with func- 
tions $ M in the integrands) is performed using an auxiliary parameter a (1 <C (j « A* ) 
(see a more detailed discussion in the next Section p| ). Detailed analysis demonstrates that 
only the region of small momenta (in comparison with the inverse electron-muon mass ra- 
tio) generates corrections of the second order in (m/M) 2 . All such corrections are connected 
with the product of the function <3>q an d of the term with factor b p in the numerator in Eq. 
(^) Q and we obtain 

2 f 1 dv / v 2 \ r 1 dz 9 , 9 3 . N /m\ 2 
S Ewv = a( Za)E F - 1 — 2 «>(l - — ) I -15^ = - a(Za) (-) E r . (7) 

Besides the electron vacuum polarization one also has to consider the muon vacuum 
polarization insertion. The respective contribution differs from the one for the electron only 
by the substitution k — * (M/m)k and has the form 




4/x^A; 4 



2 1.2 



da 2 [./n 



dk k 



+ 



k 2 + a 2 (k 2 



,2A3 



X A11 statements above are proved in Appendix [B[ 
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5E. vp = a(Za)E F ^ Cdvv 2 ^-—) f d6 sin 2 9 (2 + cos 2 6) (8) 

7T JO V 3 ' JO 



oo 1 1 p 

lo k 2 + 4 cos 2 9 4 + k 2 (l - t> 2 ) A; 2 + 16/x 2 cos 2 

Representing the last factor in the integrand in the form 

k 2 _ 16fi 2 cos 2 8 

k 2 + 16/i 2 cos 2 6 k 2 + 16fi 2 cos 2 

we easily obtain 



(9) 



7r 2 10\ a(Zce) ( m \ _ , . 



and, hence, the muon vacuum polarization does not generate contributions of order 
(m/M) 2 . Hadron vacuum polarization insertions lead to the corrections of the same order 
as the muon vacuum polarization and will be ignored in our discussion. 

III. ELECTRON-LINE CONTRIBUTION 

Radiative corrections to the hyperfme splitting induced by the radiative insertions in the 
electron line are given by the expression |19[| 

rpe-line = ^ F (_ J\ f ffi 1 ( 1 ( U] 

7T l Qfl J J m 1 (k 2 +t0) 2 U 2 + + ZO 1 ^ 



+ F-A+^o ) (74r)M v> 

where the electron factor L^ v is a sum of four terms L^ v = J2i L^ u [|19| , which after some 
transformations may be written in the form 

m - £ »*£ ( _ tf+2 t„ +aT { ^ + * kl }■ < i2 > 

Ltl = ^<7M^>(.)y ( _ fc2 + 2&fcQ + fl2)2 { c 3 A: 2 + 2 c 4 fc }, 

^ = £<W}( e) y^f ( _ P + 2 t + a2)2 { c 5 ^ 2 + 2 c 6 fc 2 fc }, 
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Auxiliary functions of the Feynman parameters a(x, y) and b(x, y) are defined by the 
relationships 



x 



y(i -y) 



1 — X 



(13) 



and explicit expressions for the coefficient functions q are collected in Table |. 
Performing contractions of the products of the electron and muon spinor structures 
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3k 2 + 2k 2 



(14) 



(7^)(e)(7^7%) = 8fc o, 
we arrive at the expression for the electron line contribution to the hyperfme splitting 



SE 



c-linc 



a(Za)Ep 



iix 2 ^ Jo dX 'o dy J m 2 \k 2 + i0) 2 



d 4 k 



1 



(15) 



1 



1 



k 2 + (i-T-ko + i0 k 2 - yr x kQ + iO 



3 /uq 2k 



cik 2 + c 2 {k 



2\2 



-k 2 + 2bk Q + a 



2 \3 



+ 



c 3 k 2 + c 4 2/c 



A; 2 + 2bk + a 2 ) 2 



3kr 



c 5 k 2 + c e k 2 2k 
■ k 2 + 2bk + a 2 ) 2 



+ 



c 7 k 2 



k 2 + 2bk + a 2 



Now we will discuss briefly calculation of the contributions to the hyperfme splitting of 
order a(Za)(m/M) 2 Ep generated by the electron factor. Let us start our consideration 
with the simplest contribution corresponding to the coefficient function c 7 in Table |. After 
Wick rotation and calculation of the angle integral we arrive at the expression 



5E = a(Za)Ep ( -) / dx dy b c 7 / dk 2 

v 7r ' Jo Jo Jo 



k 2 



(k 2 + a 2 ) 2 - Aii 2 b 2 k A 



(16) 



where the auxiliary function $i(fc) is defined in the Appendix [A]. 

As in the case of the vacuum polarization discussed above one may expand the denomi- 
nator of the integrand and obtain 



SE « a(Za)E F (-^) C dx [* dy b c 7 

V 7T ' JO Jo 



(17) 



00 u2 
dk 2 ' 



(k 2 + a 



2\2 



4 /x 2 b 2 k A 
(k 2 + a 2 ) 2 
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where the second term in the braces is at least of the second order in /i 2 due to an obvious 
inequality 4 /i 2 b 2 k 4 /(k 2 + a 2 ) 2 < 4/i 2 . 

It is easy to see that the c-integral (integral which contains function <J>f in the inte- 
grand) does not generate corrections of the second order in /i 2 since the function $f is 
/^-independent. Only the /i-integral 

5E? = a(Za)E F (— |) / dx dy b c 7 dk 2 — „M {k), (18) 
v n z/ Jo Jo Jo [k z + a z ) z 

generates corrections of the second order in the mass ratio. 

As we already mentioned in Section |T| the integrals of this kind are usually calculated 
with the help of an auxiliary parameter a, which satisfies the inequality 1 < a < yT 1 . The 
parameter a is used to separate the momentum integration into two regions, a region of 
small momenta < k < er, and a region of large momenta a < k < oo. In the region of 
small momenta one uses for simplification of the integrand the condition /ifc <C 1, and in 
the region of large momenta the condition k ^> 1. Note that for k ~ a both conditions on 
the integration momenta are valid simultaneously, so in the sum of the low-momenta and 
high-momenta integrals all cr-dependent terms cancel and one obtains a cr-independent result 
for the total momentum integral (see more detailed discussion of this method, e.g., in [f2~0|j). 
Sometimes one needs to use more stringent restrictions on the auxiliary parameter a in order 
to obtain a result with higher accuracy in /a. This well as the sample calculation 

of the /i-integral for the SE 7 contribution are discussed in more detail in Appendix |B|. The 
analysis performed in this Appendix leads to the conclusion that all corrections of order 
(m/M) 2 are generated by the low-momentum part of the /t-integral in Eq. (|i~8D , and in 
order to obtain these corrections one simply has to extract that contribution of order /t 2 to 
this integral which remains finite when the auxiliary parameter a goes to infinity. Then it 
is easy to obtain the respective contribution 

5E 7 = a(Za)E F ( - ^ ) £ dx £ dy b c 7 a = ~a(Za) {^f E F . (19) 

Calculation of the contribution to the hyperfine splitting induced by the term with the 
coefficient function C5 in Eq. (|TJ) proceeds essentially along the same lines as in the case 
of the 5E 7 contribution. Again the c-integral does not generate corrections of order /a 2 to 
the hyperfine splitting. An integral representation for the /i-integral contribution may be 



obtained from the representation in Eq. ([18]) by the substitution c 7 — > —c 5 (d/(da 2 )) 

SEt; = a(Za)E F (-^f) £ dx £ dy b c 5 jH dk 2 ^ ^(k). (20) 



Correction of order /t 2 to the hyperfine splitting is generated by the term linear in /t in the 
low momentum expansion of the function $1 (fc). Due to a higher power of the denominator 
in the integrand in Eq. (EH) than in Eq. (pT8|) the respective momentum integral is ultraviolet 



convergent and we immediately obtain 
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In the case of the weight function c% in Eq. (|15|) the //-integral has the form 

5Eg = a(Za)E F / dx dye, / dk 2 ^(k). (22) 

7T JO JO JO In* i CI ) 

Again extracting the term linear in /t from the integrand we obtain 

5£ 6 = a(Zct)£ F — y da; jf dy c 6 a = ~^ a ( Za )\J/[) E f- ( 23 ) 

Unlike all cases discussed above the /i-integral for the weight function c 4 depends not 
only on the function $1 (fc) but also on the function $2 (&) 



4^ /■! ra; roo ^2 

^ = « ^ I dx I dy b c 4 ^ cifc 2 



2$ 1 l (£;) + <$>%(k) 



(24) 



However, the conclusion of Appendix [FJ that only the low-momentum part of the /i- 
integrals leads to corrections of order /i 2 to hyperfine splitting remains valid, even in this 
case. We may repeat the considerations of Appendix [B| for this case, and check that the 
additional high-momentum contributions to the integral generated by the auxiliary high 
momentum integrals V202 and V212 (see definitions of these integrals and discussion of their 
properties in Appendix |C]) differ from the respective contributions V201 and V 2 n connected 
with the function &i(k) only by the factor 1/4. Hence they do not lead to corrections of 
order /1 2 to the hyperfine splitting. Moreover, terms linear in /i are missing in the small 
momentum expansion of the function ^ 2 {^)- Hence, one may completely ignore presence of 
the function § 2 {k) i n the integrand in Eq. (|24|) during calculation of the contribution to the 
hyperfine splitting of relative order /1 2 . Then calculation of the contribution SE^ may be 
done exactly in the same way as calculation of the contribution 5E§ above and one obtains 



5E A = a(Za)E F (-—) f dx f dy 

V IX ' JO Jo 



a 16 V ! \MJ 



(25) 



We will combine consideration of the contributions to hyperfine corresponding to the 
weight factors c\, C2, and C3. In this case, after the angular integration, we arrive at the 
expression 



SE = a(Za)Ep 



/ dx dy dk 2 k 2 < 
An 2 Jo Jo Jo 



(ci + c 2 k 2 



d \2 d 
+ 2c 3 - 



da 



'da 2 



(26) 



k 2 + a 2 



(k 2 + a 2 ) 2 - An 2 b 2 k A 



2%{k) + ^{k) 







k 2 + a 2 



da 2 > (k 2 + a 2 ) 2 - A^ 2 b 2 k A 



2$ 1 (£-) + $ 2 (£0 



Differentiation and expansion in a series in /i 2 leads to the integral 



SE « a(Za)E F 



-^r / dx I dy dk 2 k 2 < 
2n 2 Jo Jo Jo 



c 3 



Ci + c 2 k 2 

(k 2 + a 2 Y ~ (k 2 + a 2 y 



(27) 



2%(k) + $i(fc) +I2 f i 2 b 2 k i 



2ci + 2c 2 fc 2 
{k 2 + a 2 ) 5 



c 3 



(k 2 + a 



2\i 



2%(k) + ^(k) 
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-Cl 



2A/i 2 b 2 k 4 



+ 



(A; 2 + a 2 ) 3 (k 2 + a 



2\5 



2$i(fc) + $ 2 (fc) 



As usual we write the functions $i(/v) as sums of the functions <3>f (k) and $f (fc) (see 
Appendix |B|) . Let us mention here that the classical nonrecoil electron-factor contribution 
to the hyperfme splitting is produced by the term with $q = 2/(/ifc) in the integral 

above 

5E KP = a(Za)E F — ( dx j dy [ dk k 2 r 1 



(k 2 + a 



2\3 



(k 2 + a 



2\2 



(28) 



(Za)E F ( In 2 - — ). 



a 



Contributions of relative order fi 2 are connected only with the functions $f, and are 
hidden in the integral 

ci + c 2 k 2 



8E? 23 « a(Za)E F 



2i„2. 



-At / dx j dy ( dk A k 
2n z Jo Jo Jo 



2^{k) + &{(k) 



(k 2 + a 2 ) 



(29) 



^3 



(k 2 + a 



2\2 



+ 24//W$£(A;) 



2ci + 2c 2 k 2 



c 3 



(k 2 



,2\5 



(k 2 + a 



2\i 



2^(k) + mk) 



Cl 



(A; 2 + a 



2\3 



As usual, all contributions of relative order fi 2 are generated by the low-momentum 
integration region < k < cr, where one can use the small momentum approximation for 
the functions $f 



2^{k) + 



3 
2 



(30) 



2^{k) + $£(Jfe) 
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2/xA; + -(l^k) 2 . 



Due to absence of terms linear in \i in the low-momentum expansion of the function 
2$q (&) + $!(&;), the respective term in the integrand in Eq. (|29"D does not generate corrections 
of relative order /i 2 . Hence all such corrections in Eq. (p9|) are connected with the functions 
$q(/c) = and 2<&i(k) + $2^) ~ — 2fik, and are given by the expression 

2fi 2 r x f°° f C\k 4 

5E U3 = a(Za)E F dx dy dk \ lu2 ■ + (31) 



(k 2 + a 



\2b l k 



27,6 



2ci + 2c 2 k 2 
(k 2 + a 2 f 



(k 2 + a 



2\4 



-6ln2 -—)a(Za)(—) 2 E F . 
32 J v ! \M) 



Now we are in a position to write down the total contribution to the hyperfme splitting 
of relative order /1 2 generated by the electron factor, which is given by the sum of the 
contributions in Eq. (|lj), Eq. (|21]), Eq. (|23|), Eq. ([§), and Eq. (0) (see also Table [TJ 

5E el = (-61n2 - I )a(Za)(|£) V (32) 



9 



IV. MUON-LINE CONTRIBUTION 



After an obvious substitution a — * Z 2 a the electron-line factor in Eq. (O) properly 
describes radiative corrections to the muon line as well. Then the total expression for the 
muon-line contribution to hyperfme splitting of order [Z 2 a){Za)Ep may be easily obtained 
from Eq. (|ll|) and has the form 



-line 



(Z 2 a)(Za)E F 



I PX 

, , dx dy 
2n 2 Jo Jo y 



d 4 k 



m 2 (k 2 + i0) 2 \k 2 + 4/jk + iO 



1 



(33) 



+ 



k 2 - + iO 



{ ( 3k 2 - 2k 2 



cik 2 + c 2 {k 



2\2 



k 2 + 2bk + a 2 ) 3 



+ 



c 3 k 2 + c^2k 



- k 2 + 2bk + a 2 )2 



- 3fcr 



c 5 k 2 + c e k 2 2k 
( - k 2 + 2bk + a 2 ) 2 



+ 



c 7 k 2 



k 2 + 2bk + a 2 



}■ 



Apparent distinctions between Eq. ( |TT| ) and Eq. ( p3[ ) remind us that while in Eq. ( [TT]) 
the dimensionless integration momentum is measured in the units of electron mass, in Eq. 
(|33|) it is measured in the muon mass units. 

After standard transformations the integral in Eq. acquires the form 



6E n-Yme = ( Z 2 a)(Za)E F 4 F d9 sin 2 9 H dk 2 \ 1 

7T d JO JO \ 



16yU 2 cos 2 ^ 
k 2 + 16[i 2 cos 2 9 



\I(k), (34) 



where the first factor in the integrand is simplified with the help of the identity in Eq. 
1), and the second factor is given by the integral 



I(k) 



1 rx 

dx l 

o Jo 



dy{ (2 + 



cos 



(ci sin 2 9 + c 2 k 2 



d \2 



da 2 



+ 2c, 



d 

da 2 



(k 2 + a 2 ) (35) 



d d 
- 8bcJ2 + cos 2 9)cos 2 9— + 12&cos 2 flf c 5 — - c 7 ) 

an 2 V rtd 2 / 



12c 6 cos 2 ^ ( , 2 + a2) }____^ 



4& 2 A; 2 cos 2 



The first factor in the parenthesis in the integrand in Eq. fl3~j| ) generates only corrections 
of relative order /i, since the integral I(k) is /i-independent. The respective correction 
to hyperfme splitting was calculated long time ago |T7tp5lp6| , and we will not consider it 
here. Correction of order fi 2 in Eq. ([34]) could be connected only with second term in the 
parenthesis which apparently generates a contribution of order /i 3 . However, the integral 
I(k) behaves like 1/k at small momenta, and this singular behavior reduces the apparent 
power of /i, generating a correction of relative order /i 2 . In order to obtain this correction 
we have to calculate the leading term in the low- momentum expansion of the integral I{k). 
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Let us outline briefly the respective calculation. First, we note that the most singular 
k behavior of the integral in Eq. (|35|) is connected with the region of small values of the 
Feynman parameter x. Hence, in the limit of small k and x we preserve in the integrand 
only nonvanishing terms. The calculation is facilitated by the change of variable y = xz. 
Then use the following approximations 



a 2 « 



x 



b « 1, k 2 + a 2 



x 



(36) 



c x ss 16 ( - - 3 -21nz), c 3 



1 , 1 



« — ( 26 + 6z — 16 In x) 

x z 



c 4 ~ 5, c 5 

z 



8, c 6 



l-^ 2 
, c 7 ^ - 



X 



X 



The term with the coefficient function C2 in Eq. (|3^) may be omitted since it is explicitly 
multiplied by the small factor k 2 . Then the leading infrared contribution to the function 
I(k) may be written as 

J(fc) « / d^ 2 / dxf (2 + cos 2 fl)[ 16sin 2 ^( --3-2\nx)x(— ) x (37) 
JO Jo 1 z K OX J 



+ 2(- -26 + 6^- 16 \nx)-0-x } + 8(2 + cos 2 9) cos 2 0(-2 + 5z)x-^- 
z ox 1 ox 



+ 12 cos 2 9 



s 9 

(Q-8z)x- 2 

Ox 



d 



+ 12 cos 2 0(1 - z)— x } 

(IT, J 



dx > x 2 + Ak 2 z 2 cos 2 
Integrating by parts we reduce all infrared-hard terms to the standard integral 



dx 



IT 



lo x 2 + Ak 2 z 2 cos 2 9 4fczcos0' 
and arrive at the relationship 



(38) 



/(*) 



dz z A 



32(2 + cos 2 0)(1 -zsm 2 9) + 8(2 + cos 2 9) cos 2 9(2 - 5z) (39) 



12 cos 2 0(6 -8z + 2) 



7T 



4tt 2 + cos 4 



cos 3 A; cos 9 



Now we easily obtain radiative-recoil correction of the relative order fi 2 from the second 
term in the integrand in Eq. ([34]) 



6E 



ml 



17 
12 



M) 



Er 



(40) 
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V. DISCUSSION OF RESULTS 



Let us consider first the electron-line contribution in Eq. This result demonstrates 
that the discrepancy between the results in |17j and in [jl8| , |l9| for the radiative-recoil cor- 
rections induced by the electron line cannot be explained by the contributions of the terms 
of higher order in the mass ratio to the result in ||17|| . Numerically we obtain from Eq. (|3~2"D 



5E el « -0.03kHz, (41) 
which is too small to explain the difference 0.22 kHz between the results in |T7] and 



T8|]19|| . We would like to remind the reader that the coinciding results in |18| and [19| were 



obtained completely independently in different gauges, so there is little, if any, doubt that 
this result is correct. 

The total one-loop radiative-recoil correction of relative order (m/M) 2 is given by the 
sum of the results in Eq. ([?[), Eq. (|32|), and Eq. 



^-^rad-rec 

or, recalling that Z = 1 for muonium 



61n2 - ^ ) a(Za) - — (Z 2 a)(Za) 



m \ 2 
MJ 



E F , (42) 



SE Tad . iec = (-6 In 2 - — ) a 2 J Ep « -6.32555 ... c, 2 (-) Ep. (43) 

Let us mention a minor subtlety connected with the definition of the Fermi energy in 
Eq. (fi"2"|). As is well known, there are two common definitions for the Fermi energy: one 
includes the muon anomalous magnetic moment, and the other does not. Both definitions 
are useful in presenting results for the high order corrections to the hyperfme splitting. 
Corrections, generated at the scale of order 1/m or larger (small intermediate momenta) are 
multiplied by the Fermi energy which includes the muon anomalous momentum. But for 
the corrections generated at the scale comparable with the muon Compton length 1/M the 
muon anomalous momentum cannot be factorized, and such corrections are multiplied by the 
Fermi energy without the muon anomalous moment. As we have seen above, radiative-recoil 
corrections of the second order in mass ratio originate from the distances larger than the 
electron Compton length, and, hence, the Fermi energy in Eq. (|42"D should include the muon 
anomalous moment. In this way we properly take into account some corrections of higher 
order than a(Za)(m/M)Ep but due to the smallness of the correction under consideration 
this is irrelevant from the purely phenomenological point of view. 

Numerically the contribution to hyperfme splitting obtained in this paper is equal to 

^rad-rec = -0.0351... kHz, (44) 

and should be taken into account, along with some other small corrections to hyperfme 
splitting calculated recently, when comparing the pending experimental results with the 
theory. 
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APPENDIX A: STANDARD AUXILIARY FUNCTIONS 



All contributions to hyperfine splitting in the main body of this paper are written in 
terms of the auxiliary functions (n = 0, 1, 2, 3) defined by the relationship 



l r 



TT/jl* JO 



d9 sm 2 9cos 2n 9 



(k 2 + a 2 ) 2 - 4fi 2 b 2 k 4 



( k 2 + /i- 2 cos 2 9 ) [ (k 2 + a 2 ) 2 + 4 b 2 k 2 cos 2 9 



(Al) 



The integral over angles may be explicitly calculated, and the result of the integration is 
conveniently written as a sum 



where 



(A2) 



(A3) 



(A4) 



+ 16' 



2' 



(A5) 

(A6) 

(A7) 
(A8) 
(A9) 



(A10) 



*?(*) = ic[ic(icw{ic) - \) + 



1 1 



1 

16' 



(All) 



The standard function W(£) has the form 
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W(0 = - 1. (A12) 

and 

«, = >=\ ic = (A13) 
One may easily obtain asymptotic expressions for the function W(£) 

(A14) 

lim - 

High- and low- momentum asymptotic expressions for the functions <&i(k) also may be 
easily calculated. Let us cite low-momentum expansions, which where used in the main text 
for calculation of the contributions to the hyperfme splitting of relative order fi 2 

« " 1 + Y' (A15) 

&t(k) « ^ - M + (/1A;) 2 , (A16) 

m) ( A17 ) 



APPENDIX B: CALCULATION OF ^-INTEGRALS 

In this Appendix we will consider the main tricks used in calculation of the //-integrals. 
Consider as an example calculation of the integral relevant for the 8E7 correction 

rl j-x roo 

m= J,*/, *»><>, J, dejp+^w). (bi) 

We need to calculate contributions order \i to this integral. As was discussed in the 
main text we introduce an auxiliary parameter a (1 C <r < and consider separately 
contributions of the large and small integration momenta 

F(jm)=F < (jj.) + F > (jm). (B2) 

In the large momenta region we would like to use the relationship k > a ^> 1 in order to 
simplify the integrand before integration over momenta. The simplest approach would be to 
expand the denominator in a power series over a 2 /k 2 prior to integration over the Feynman 
parameters x and y. However, it is easy to see that the second term of this naive expansion 
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would lead to a divergence in the integration over y when y goes to zero. Of course, the 
initial integral is convergent, and this divergence is a result of an improper expansion. A safe 
approach would be to calculate first the integrals over the Feynman parameters, and only 
then to expand in large k 2 . Emergence of the logarithm of k 2 in this expansion corresponds 
to the divergence which was encountered in attempting to make the large k 2 expansion prior 
to the Feynman parameter integration. 

However, we would like to avoid tedious calculation of all integrals over x and y prior to 
the large momentum expansion. As we will demonstrate now this may be achieved with the 
help of a trick, which will significantly simplify further calculations. The idea is to use the 
following representation for the coefficient function in Eq. flB~T|) 



bc 7 = bc 7 + {bc 7 — bc 7 ), (B3) 

where 

ter=(l- 2y)(b(x,y)c 7 (x,y))\ y=0 . (B4) 

The weight function (bc 7 — bc 7 ) is obviously proportional to y, and, hence, the high- 
momenta integrals with such weight function admit naive expansion over 1/k 2 prior to 
integration over the Feynman parameters. Integration over y in the integral with the weight 
function bc 7 is simplified due to the presence of the factor 1 — 2y = [y(l — 2y)]'. 

Let us consider first the integral with the weight function bc 7 

F(jx) = 2 dx dy { - L(i _ 2y ) / dk 2 (B5) 



o Jo x Jo (k 2 + a 2 ) 2 

In order to calculate the high-momenta contribution we first integrate over variable y, 
expand the result in the series over 1/k 2 , and then integrate over x 

^ = 2 /> 2 {i^- 3^}*«*>' < B6 > 

Note that we took into account only the first two terms in the expansion over 1/k 2 . As 
we will see below, due to the special choice of the auxiliary parameter a, these terms are 
sufficient for calculation of the integral above with accuracy linear in /x. The remaining 
momentum integration is easily performed with the help of the standard integrals Vkmn 
which are calculated in Appendix 0, and we obtain 

F > fr) = 2(\v wl -\v B1 ) «2[ - j In (2/i) - ~ - \\na (B7) 



+ — - — -lna H 

2 4 a 2 V 6 12 > J 

The result above is a power series in the small parameters (ficr) and l/o 2 . In order to 
get rid of higher terms in the expansion we further specify the magnitude of the auxiliary 
parameter a ~ 1/ JJi. With such a choice of this parameter we need only explicitly written 
terms in the relationship above in order to obtain the value of the integrand with linear 
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accuracy in /x. The same specification of this parameter also justifies consideration of only 
the two first terms in the expansion over 1/k 2 above. 

Calculating the low momentum integral we first expand function <&i(k) in the small 
parameter \ik (up to and including the term (/ifc) 2 < (/^cr) 2 ) 

&£{k) « 5 - ji* + (jA) 2 + 0((^) 3 ), (B8) 
and then perform the momentum integration 



era a 

2cr + — — 3a arctan — 

a 1 + er a 



2 ° 2 a 2 2l ° 2 + a 2 
a — — — 2a in 



a 2 + a 2 a 2 



Calculating remaining integrals and omitting the higher order terms with the help of the 
condition a ~ 1/ JJi we obtain 



1 1 , 15vr 2 uo ua 2 1/1 1 

1 lncx H a - - h ^—^ 1 -lntr H 

48 4 128 P 2 4 a 2 \ 6 12 



(BIO) 



All terms with the auxiliary parameter a cancel in the sum of the high- and low- 
momentum integrals and we obtain with linear accuracy in the small parameter \i 



1 , ,~ x 7 



2 v r/ 24 



15tt 



2 



F(n) « - -In (2^) - — + — — /x. (Bll) 
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The first term in this relationship gives the contribution of the first order in the small 
mass ratio to hyperfine splitting, and is of no interest to us. The second term leads to the 
correction of order (m/M) 2 and calculation of such terms is the main goal of this work. 

Let us consider further the origin of this term. First, it is not connected with the high- 
momentum integral, which generated only corrections of the previous order in the mass 
ratio. The term which generates the correction under consideration in the low-momentum 
integral in Eq. ( |B9|) and Eq. ( piO|) may be easily identified. This is the term which explicitly 
contains factor /i generated by the low- momentum expansion (see Eq. (|B8j) ) of the function 
§i(k), and moreover it is the only term of order fi in the integral which remains finite when 
a goes to infinity. We may be even more concrete. It is easy to see that the only term which 
remains finite when the upper limit in the integral Jq dk 2 k 2 /(k 2 + a 2 ) n /x/c goes to infinity 
has the characteristic form of arctan(<r/a). 

Now we are ready to describe a recipe for obtaining the term of order \x in the integral in 
Eq. (|B5| ) almost without calculations. First, we have to calculate only the low- momentum 
contribution to the integral in Eq. ([B~5|), and consider in the integrand only the term linear in 
[i generated by the respective term in the low- momentum expansion of the function 3>i(fc). 
At the second step we calculate the momentum integral and preserve in the result only 
contribution of the arctangent which gives factor tt/2. At the last step we have to calculate 
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the remaining double integral over the Feynman parameters. It is easy to check that this 
recipe immediately reproduces the term linear in /x in Eq. (|B11| ). 

Up to now we have calculated only the linear in /x contribution to the integral F{\x) in 
Eq. (p5|) . However, it is easy to see that the recipe for extracting linear in /x contributions 
which we have just formulated is equally applicable for the total integral F(fi). The only 
difference between the integrals with the weight functions bc-j and bc 7 is in the details of the 
calculation of the high momentum contribution to the integral. The representation of the 
weight function in the form of the sum in Eq. (|B3|) , simplified the calculation of the high- 
momentum contribution generated by the weight function bc 7 — bc 7 , but as we have just seen 
the high-momentum contribution does not generate linear in /x contributions to the integral 
F(n). Hence, we may apply the recipe above for calculation of the total contribution of 
order /x to the integral F(/j,), the only difference being that the role of the weight function 
will play now the total weight function bc 7 . 

It is clear that the same recipe will work each time when the /x-integral for the contribu- 
tion to the hyperfine splitting contains function $>x(k) in the integrand. Moreover, one may 
check explicitly that the recipe is valid even when the /x-integral depends on other functions 
$o 2 (A;), as happens for some contributions generated by the electron factor. All respective 
calculations in the main body of the paper are made with the help of this recipe. 



APPENDIX C: STANDARD HIGH-MOMENTUM INTEGRALS 

As was mentioned in the main text and Appendix |B| all high-momentum contributions 
to hyperfine splitting may written in terms of the standard integrals 

r°° dk 2 

Vlmn = I J^y^r KM, (CI) 

where I — 1, 2, m — 0, 1 and n = 1,2. 

These integrals should be calculated with linear accuracy in /x. As was discussed in 
Appendix [B] this means that integrals should be calculated up to and including small terms 
of order (/xcr) 2 (we remind the reader that o « l/y^x). 

Consider in more detail the integrals with n = 1 



oo Ah 2 

'"' . . .. . , 2 ; 2 



(*») 



Vimi = / 77™ (In k) m - nkJl + fiW + fi 2 k 2 + - . (C2) 



First, we change the integration variable 



2 



t = ^l + fi 2 k 2 - fik. (C3) 
In terms of this new variable typical integrals acquire the form 

- 1 * 4^- < c4 > 

rU + (I + t 2 ) 1 - t 2 
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^2oi = V \ ta dt ^"ff , (C6) 





where t CT = VI + A^ " 2 ~ A* " ~ 1 — ^ + (A* 20 " 2 ) /2- 
Let us, for example, calculate the integral Vm 

= r*'(r§F- 1 )'»^r s ^ + v "- < C8 > 



The first term in Eq. (|CS|) may be written in the form 



V = - / <<T rfln(l-t 2 )ln^ — - (C9) 
Jo 2/it 



1 dt 
= ~ + In (2/^) In (1- A 2 ) - j o y ln(l-t 2 ). 

With the help of an exact relationship 

1-A 2 = 2fiat a , (CIO) 

V may be written in the form 

V = ^ln 2 (2/i) - -In 2 a + iln 2 t CT + ln(2/i)lnt CT (Cll) 

t Ha t 

Expansion over (/xcr) is facilitated by the relationship 

lnt a « - //o- + 0[ (/iof ], (C12) 

and we obtain 

V « ^ln 2 (2/i) - iln 2 a + g + (^)( ln<r - 1 ) + 0[ (/ia) 3 ]. (C13) 
Approximate expression for the integral V" may be obtained in the same way 

V" « i + ^ln(2/x) + (H( Ino- - 1 ) + (//of( - \na + \ ), (C14) 
and then we obtain an approximation for the integral in Eq. (|C8f) 

Vm ^ ^ln 2 (2/i) + iln(2/i) + ^ + I - IhiV (C15) 

+ (//<r)( 21n<r - 2) + (/i<r) 2 ( - lno- + - 
Approximate expressions for other integrals Vi mn may be obtained in the same way, and 



they are collected in Table |T] . 
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TABLES 

TABLE I. Coefficients in the Electron-Line Factor 



Cl 


16 


(1 — x)(x — 3y) — 2ylnx 


2/(l-2/) 3 


C2 


4 


- (1 - x)(x - y - 2y 2 /x) + 2(x - 4y + 4y 2 /x) lnx 


y(^-y) :s 


C3 


1 


1 - 6x - 2x 2 - (y/x)(26 - 6y/x - 37x - 2x 2 + 12xy + 16 lnx) 


y{i-y)' 2 


C 4 


i 


2x - 4x 2 - 5y + 7xy ) 
V / 


2A 1 y) 


C5 


i 


' 6x - 3x 2 - 8y + 2xy ) 




C6 


ft2 x-y 




C7 


2^ 
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TABLE II. Electron-Factor Contributions to Hyperfine Splitting 



SEi 


a(Za){^) 2 E F 




(-6 In 2-|) 


5E 4 


3 

16 


SE 5 


3 

32 


5E 6 


3 

32 


5E 7 


9 

32 




(-6 In 2 - |) 
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TABLE III. Approximate expressions for the integrals Vi mn 



^100 


21n(2/x) - 2 + 2W - (fia) 




V 110 


- ln 2 (2^) + 21n(2/x) - ^ - 2+ ln 2 <7 + 


-In (j + 1 ) 


Vioi 


- ln(2/x) - i - W + 2(/x<r) - (^a) 2 






\ In 2 (2/i) + \ In (2/i) + fJ + ±- ±ln 2 a + (/xa) (2 In a - 


2) +( /[ ia) 2 (-lna + i) 


V20O 


1 

~^ 




V201 


1 

2^ 




V202 


1 

8^ 




V210 


-£■( In a + \ ) 




V211 


2M In a + 




V212 
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